New Constraints on the Highest-Energy Cosmic-Ray Electrons and Positrons 
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At energies above a few TeV, no measurements of the cosmic-ray electron spectrum exist yet. By 
considering the similarity of air showers induced by electrons and gamma rays as seen by ground- 
based arrays, we use published limits on isotropic gamma-ray fluxes to place first constraints on the 
> 10 TeV electron spectrum. We demonstrate that, due the proximity of known sources, the flux of 
such electrons (and positrons) can be large. We show how these smoothly connect to lower-energy 
positrons measured by PAMELA and relate to exciting new indications from Fermi. 
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Introduction. — Interest in the cosmic-ray electron 
spectrum at Earth is at an all-time high, arising from 
both astropbysical and dark matter-related concerns [l|, 

0, 1 I, !, 1 0, S S E3- In the GeV ener §y ran S e , 
much greater clarity than what had existed in preced- 
ing decades has been brought by the PAMELA [l| and 
Fermi 2] space missions. 

Such measurements, which observe electrons and 
positrons directly, become difficult above ~ 1 TeV due 
to their fixed detector areas relative to declining particle 
fluxes. Using the indirect technique of observing atmo- 
spheric air showers, HESS has pushed the energy fron- 
tier up to several TeV [3|, although above this, no elec- 
tron measurements have been reported. With the histor- 
ical lack of local multi-TeV gamma-ray sources and the 
soft spectra of secondary e ± resulting from p-p scatter- 
ing , little if any signal may have been expected. 

Due to the similarity of the electromagnetic showers 
produced in the atmosphere by energetic electrons and 
gamma rays, they are nearly inseparable via ground- 
based observations [H],[l3j]. Here, we translate published 
limits on isotropic gamma-ray fluxes (e.g., [lJ HH HE 
[lH) to constrain the e~ + e + spectrum to > PeV ener- 
gies. 

We further show that our derived limits are relevant in 
light of recent measurements by Fermi and TeV gamma- 
ray telescopes, particularly the discovery and detailed ob- 
servations of high-energy sources both near and far that 
indicate e^ production in the GeV and TeV regimes. 
Considering these observations, with an improved ana- 
lytical treatment of e ± propagation to calculate expected 
source contributions at Earth, we conclude that it would 
be surprising if the influence of a nearby pulsar was not 
present in cosmic-ray positron data. 

New Cosmic-Ray e^ Limits. — While the ~ 1 m 2 Large 
Area Telescope (LAT) of Fermi [2j has brought a sharper 
picture of the e" + e + spectrum up to 1 TeV, to progress 
further requires a much larger effective area for particle 
collection. This can be accomplished by, in lieu of direct 
particle identification, examining the showers that result 
when energetic particles scatter in the upper atmosphere. 

To observe an electromagnetic shower (initiated by a 
gamma ray or electron) requires rejecting the large back- 



ground due to cosmic-ray protons. Ground-based air 
Cherenkov telescopes (ACTs) operate by imaging the 
shower and observing the differences between hadronic 
and electromagnetic cascade development [17j . which al- 
lowed HESS to measure e~ + e+ up to several TeV 
(Fig.[T]). Alternatively, detectors that operate by directly 
detecting the long-lived products of air showers as they 
reach the ground make use of the low muon content of 
electromagnetic showers relative to hadronic events (l7j . 

Searches for isotropic fluxes of gamma rays have been 
conducted over a wide ran ge o f energies by examining 
muon-poor showers fbH . HH Il6|] . While these have typi- 
cally only resulted in upper limits, they can be valuable. 
Since a gamma ray first produces an initial e pair to 
begin a cascade, its shower will look very similar to that 
of an electron of equivalent energy, only located ~ one 
radiation length deeper in the atmosphere [12]. HESS 
electron data can thus be regarded as upper limits on an 
isotropic TeV gamma-ray flux [E [l8j], while air shower 
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FIG. 1: The cosmic-ray electron spectrum at Earth. Shown 
are direct e~ + e + measurements from Fermi [2L measure- 
ments of e~ + e + based on showers from HESS 0, 0], and a 
baseline oc E~ 3 . These can be compared to our limits derived 
from gamma-ray experiments [13 . Il4l . [l5l . [l6l | at > 10 TeV. 
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arrays do not have hope of exploiting this subtle differ- 
ence, rendering diffuse gamma rays and electrons insep- 
arable (l3| . 

By using isotropic gamma-ray limits from air shower 
arrays, we derive new limits on the cosmic-ray elec- 
tron spectrum in a regime currently lacking constraints. 
These limits are often quoted as the fraction of mea- 
sured gamma-ray to proton intensity, J 7 //cr- We use the 
cosmic-ray nuclei spectrum [19( in Fig. [TJ In the range 
10- 10 5 TeV, we obtain electron limits from HEGRA pi, 
CASA-MIA Q, GRAPES [H, and KASCADE [r| 
data, leading to the limits in Fig. [TJ 

We have conservatively assumed that the electron 
fraction of the electromagnetic showers seen is 100%. 
It is likely that dedicated re-analyses of the data can 
strengthen these constraints, and an analysis of Mila- 
gro [20] data could probe the 10-100 TeV region. Even 
so, we see that our limits are already competitive when 
compared to previous measurements at lower energies. 

Nearby e Factories. — Due to their short lifetime 
against radiative losses, any electron flux measured at 
E > 10 TeV must be produced in the very-recent his- 
tory of a nearby source. For this reason, we will con- 
sider sources within ~ 500 pc that still exhibit evidence 
of particle acceleration to very-high energies, using more 
distant sources for added guidance. Observations, prin- 
cipally of extended gamma-ray emission, have recently 
revealed good reasons to believe that e ± are being pro- 
duced, accelerated up to multi-TeV energies, and escap- 
ing from the high-energy-density environments of pulsars: 

• HESS has intensively observed the distant (~ 4 kpc) 
pulsar wind nebula (PWN) HESS J1825-137 ;21]. This 
included untypically-long exposure times due to the 
PWN being serendipitously within the field-of-view of the 
variable microquasar LS 5039, for which repeated mea- 
surements of the light curve were taken. These revealed 
an extended wind of > 10 TeV e ± containing > 10 48 erg 
reaching > 100 pc in only ~ 20,000 yr [2l[. This indica- 
tion of a pulsar "TeV mode" carrying > 10 48 erg of e 
to great distances likely would have been missed without 
these deep observations, and thus may be common. 

• The proximity (290 pc) and relative youth (~ 
11,000 yr) of the Vela pulsar has permitted detailed ob- 
servations of its PWN, Vela X. There, HESS inferred a 
population of e ± with energies reaching ~ 100 TeV (22j. 
Modeling source e ± injection spectra as dN/dE cx 
E~ a e~ E ' Ema * , the HESS measurements imply a = 2 and 
^max = 70 TcV for Vela X [22j • Intriguing evidence 
has also been obtained at lower energies, with a distinct 
population of multi-GeV, radio-synchrotron emitting 
(with an E~ ls spectrum) modeled in Ref. 23]. The 
(preliminary) Fermi discovery of extended GeV gamma 
rays in Vela X confirms this and, importantly, indicates 
a high-energy cutoff in the spectrum at ~ 130 GeV J24[. 
The amplitude of this signal implies that a pulsar "Ge V 
mode" can also produce > 10 48 erg of e . This suggests 
that the total composition of Vela X includes ~ 100 GeV 
e ± accelerated by the pulsar itself — likely associated 



with pulsed GeV gamma rays — and TeV e ± from shock 
acceleration of the pulsar wind. As observed, the GeV 
component contains ~ 100 times more energy, although 
multi-TeV particles may have already exited the system. 

• The discovery of extended ~ 35 TeV gamma-ray 
emission by Milagro [2C| surrounding the nearby (~ 
200 pc) pulsar Geminga indicates a close, active source 
of e [5j. These data imply a firm lower limit on the 
maximum particle energy of > 100 TeV, and may ap- 
proach 1000 TeV. In all of these systems, the implied par- 
ticle multiplicities needed to account for the gamma rays 
via inverse-Compton scattering require e pair produc- 
tion 0, [25[. Since, like Vela, Geminga possesses bright 
pulsed GeV emission, it may as well have resulted in an 
abundance of ~ 100 GeV e ± . While Vela is too young for 
GeV e ± to have reached us yet (as we will soon see), the 
greater age of Geminga (~ 300,000 yr) may allow for a di- 
rect test of the commonality of dual high/low-energy pul- 
sar e populations. To determine this, we first address 
the propagation of highly-energetic e ± in the Galaxy. 

TeV Propagation. — In spherically symmetric ge- 
ometry, the diffusion equation governing the particle den- 
sity at a given location/time/energy, n(r, t, E), is [2III27I] 



On D(E) 8 2 dn 3 

m = —d- r r Tr + dE [h{E)n] 



Q 



(1) 



with energy losses parametrized as b(E) = —dE/dt, dif- 
fusion coefficient D(E), and source term Q. We first 
consider a single burst from a point source with a gen- 
erated particle spectrum dN/dE g , so that Q(r,t,E g ) = 
5(r) 5(t) dN/dEg. Using the Syrovatskii propagator [261 ] 
(as in [Elll ) yields 



n s {r,t,E) = 



-r 2 /r 2 di , dN dE 



* 3/2 r% f dE g 



dE 



(2) 



where dE g /dE maps the energy at generation E g (E, t) to 
the observed E after losses. In the high-energy regime of 
interest here, most relevant are inverse-Compton losses 
on the CMB (energy density ~ 0.3 eVcrn" 3 ) and syn- 
chrotron losses due to a ~ 5 /xG magnetic field (~ 
0.2 eV cm" 3 ), so that b(E) = b E 2 with b 
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Integration yields 1 / E — 1/E g 
so that dE g /dE = (E g /E) 2 . 

The diffusion radius is rdif(E,t) = 2y/X(E, t), with 



X(E,t) 



dt'D[E(t')] 



dE 



, D(E>) 
b(E') 



(3) 



We parametrize the diffusion coefficient as D(E) = 
do(l + E/Eq) s , with Eq ~ 3 GeV, and consider values of 
d = 2 x 10 28 , 4 x 10 27 cm 2 s- 1 , 5 =0.4, 0.5 [UHil cor- 
responding to "fast" and "slow" models. Since D(E) ~ 
do(E/Eo) 5 at energies E > Eq, the diffusion radius sim- 
plifies to 



r dif (E,E g 



d E 
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TABLE I: Parameters used for the two Geminga scenarios. 
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FIG. 2: Our modeled contributions to the e + /e~ spectra. 
Bottom: Shown are e~ +e + measurements from Fermi [2] and 
HESS 0,01; direct e" from AMS 0; our fit to these; and our 
> 10 TeV limits from Fig. [T] For Geminga, we show a model 
in which the TeV mode dominates and one with comparable 
GeV/TeV modes. Top: A comparison of the expected cosmic- 
ray positron fraction from the two models with data up to 
~100 GeV 0, S2,|i3. 



rate dN/dE g , as Uq (E) = f*dt'h(r,t\E). Note that for 
sources such as pulsars, r can vary with t 

Beyond the Positron Excess. — Using the Jiittner for- 
malism, we can better estimate the extent to which a 
source may influence the spectrum at Earth. For mag- 
netic dipole braking, pulsar spin-down power evolves as oc 
(l+£/i )~ ( " +1)/( ™~ 1) , with n = 3 and a pulsar-dependent 
timescale, io- Assuming that the e^ output is simply pro- 
portional, the corresponding evolution will be C e ± (t) = 
[E G /t G ) [1 + (t G - t)/to] -\l J tG dt' [1 + (t G - f)/t ] - 2 . We 
normalize to the e^ luminosity of the source, C e ± (t) — 



J E EdN/dEdE (with E m 
For Geminga, E ~ 10 34 5 



1 GeV). 

erg s _1 at present, and we 
use to ~ 3 x 10 4 yr (corresponding to a spin period at 
birth of ~ 70 ms; see [3J] for how this can vary). We con- 
sider two distinct scenarios: a case from Ref. [5[ in which 
the TeV mode dominates and another with comparable 
energy in the dual GeV and TeV modes. For the TeV 
spectra, we assume E maK — 200 TeV, and a = 2 [5J. For 
the GeV spectrum, we use parameters consistent with 
Vela X: £ max = 150 GeV, and a = 1.8. 



We use rdif(E,E g ) to gain insight into which sources 
could contribute to the multi-TeV electron spectrum. If 
a particle generated with even E g = 100 TeV is to be 
detected at E — 10 TeV, it may propagate a maximal 
distance of r^if ~ 160 (230) pc in slow (fast) diffusion 

models. Note that the term e~ (r /r ^V r Ii/ in Ec l- © 
strongly favors particularly nearby sources. 

Closer inspection of the Syrovatskii solution reveals 
that more care is needed at these high ene rgies. For 

Eg ~ E, E ~ b t and r dlf (t) ~ 2 y / D(E)t, the quan- e +/(e+ + e~). For the denominator, we directly use a 



In Fig. [21 we display the local flux of e~ + e + , J© = 
(c/47r)n0, for these models, as described in Table [T] 
These require ~ 40% of the spin-down power be con- 
verted to high-energy pairs, within the range inferred 
from Vela X [23j . We also examine the expected ratio of 



tity Vdif{t) = rdif(t)/t ~ 2y / D(E)/t can exceed the 
speed of light. This is the known problem of "super- 
luminal" diffusion, which lies in the fact that the diffu- 
sion solution is not relativistic, similar to the Maxwell 
distribution yielding values in excess of c [13, HH HH ■ A 
phenomenological resolution was discussed in the context 
of ultrahigh-energy cosmic rays in Ref. [33j , which uses a 
propagator based on the Jiittner particle distribution (3(| 
to explicitly limit fluxes to v < c, while preserving diffu- 
sive behavior at lower energies. Now, 



nj(r,t,E) 



a 



e- a l^- 



o 



dN dE n 



ATr(ct) 3 (l-£ 2 ) 2 Ki(pt)dEg dE ' 



(5) 



where £(r, t) — r/ct, 9 is the step function, K\ is the mod- 
ified Bessel function, and a{E,t) = c 2 t 2 /(2\(E, t)) [pj. 

This formulation has the added benefit of removing 
spurious features appearing in the non-relativistic solu- 
tion for source spectra harder than E~ 2 and can easily be 
generalized to a continuously-emitting source at distance 
r and time t with a time-dependent particle injection 



fit to the measured data (neglecting the final HESS da- 
tum), as shown in Fig. [5] In the numerator, we include 
secondary fluxes [ll| to match low-energy data. As we 
see in the top panel, the two models can easily diverge 
at ~ 100 GeV. A more GeV-dominated case would drop 
more dramatically. 

A few remarks are in order: (1) The underlying Galac- 
tic e~ component is not known, but must be cut off at 
some point in order to not overshoot high-energy data 
(as seen in Fig. This must make up the difference 
between any model and the full measured e~ + e + spec- 
trum. (2) Multiple nearby sources may contribute, al- 
though Geminga and Vela are the only bright gamma- 
ray pulsars within ~ 300 pc. (3) Low energies are more 
sensitive to the source's initial position, while at high en- 
ergies the present position is most relevant (which can be 
directly measured) . (4) The emission properties (spectral 
index, cutoff, pair conversion efficiency, etc.) may evolve 
in time, while the highest-energy emission can reasonably 
be tied to what is seen today. The observations of Vela 
do give hope of examining lost history, though. 
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Conclusions. — The flux level reached by our derived 
limits is encouraging for the prospects of upcoming gen- 
erations of electron experiments, and already constrains 
the presence of a nearby, very-high-energy source. Mean- 
while, our handling of propagation in the presence of en- 
ergy losses eliminates the appearance of spurious super- 
luminal solutions. This serves as a step towards a more 
accurate description of time-dependent e ± propagation. 
We have considered the nearby pulsar Geminga, detected 
in multi-TeV gamma rays, as the best motivated source 
for multi-TeV electrons reaching Earth today, although 
Vela X may also contribute contingent on escape. 

A nearby supernova remnant could result in elec- 
trons [35! ] . although age and distance play significant 
roles, with the recent history of acceleration being quite 
important due to the short cooling time at high energy. 
Prominent SNRs include [H the Vela SNR (-250 pc), 
Cyngus Loop (~ 440 pc), Monogcm Ring (~ 300 pc), and 
Loop I (~ 100 pc). Due to their large angular extents, 
it is difficult to establish with ACTs whether these are 
active. Searches via synchrotron radiation or wide-field 
gamma-ray instruments are better suited. The only SNR 
potentially within 500 pc and detected in TeV gamma 
rays is Vela Junior, with an uncertain distance ranging 
from ~ 200 pc (corresponding to an age of only ~ 500 yr) 
to — 1 kpc (~ 5000 yr) [37j , while the "Boomerang" 
SNR/PWN at ~ 800 pc has been seen by Milagro 

If features in the e~ + e + spectrum, such as the change 
in slope at ~ 100 GeV seen by Fermi or the drop at 
— 1 TeV measured by HESS, are due to a transition be- 
tween different sources, classes, or populations, these 
should correspond to features in the positron fraction. 
Measurements of the separate e _ /e + spectra will be vi- 
tal in determining the Galactic component and in dis- 
tinguishing SNR, dark matter [381 ] . and pulsar contri- 
butions. For regions of pulsar domination, the positron 
fraction should saturate at ~ 50%, while SNRs result 
entirely in primary electrons (unless secondary accelera- 
tion occurs @]). While PAMELA can reach 300 GeV for 
e+ [li], AMS-02 will go to ~ 1 TeV 0, and ACT mea- 
surements using the varying position of the moon shadow 
for e + /e~ may reach several TeV [4l|. These combined 
observations give hope for connecting very-high-energy 
emission to lower energies to build a complete picture of 
the spectra of cosmic-ray electron and positrons. 
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